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From Navier-Stokes turbulence numerical simulations we show that for 
the extended self similarity (ESS) method it is essential to take the third 
■ order structure function taken with the modulus and called D^{r), rather than 

the standard third order structure function D^{r) itself. If done so, we find 
ESS towards scales larger than order ~ lOr/, where r] is the Kolmogorov 
bJQI scale. If D2,{r) is used, there is no ESS. We also analyze ESS within the 

Batchelor parametrization of the second and third order longitudinal structure 
function and focus on the scaling of the transversal structure function. The 
. Re-asymptotic inertial range scaling develops only beyond a Taylor- Reynolds 

\ number Rex > 500. 

> 

g ■ I. INTRODUCTION 

O , Extended self similarity (ESS, [|l]-0]) has been most useful in determining scaling ex- 

Q>^ \ ponents in experimental and numerical turbulent flow. In ESS, the p*'^ order longitudinal 
\ velocity structure function 

^ : Diir) = {[{uix + r) - uix)) • e^f) (1) 

O . 

J2 ■ is plotted against the third order structure function. Here, is the unit vector in r direction; 

the unit vector used below is perpendicular to e^. The original motivation for picking the 
third order structure function was the Howard-v. Karman-Kolmogorov structure equation 



- 
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D^{r) = -^er + 6u^D^{r), (2) 

saying that in the inertial subrange (ISR) D^{r) scales like oc r and therefore Dp oc 
r'^p oc (Dg (r))^p has the same scaling exponent as a function of r or of D^. However, because 
of the poor statistical convergence, rather than (r), the third order structure function 
D^^{r), calculated with the modulus of the velocity difference, is taken and it is argued 
that D^^{r) would also scale linearly with r in the ISR 0. The resulting exponents, 
which in general have to be distinguished from the 's are denoted as , defined by 

oc {Df^p , and they are found to be remarkable universal, i.e., independent of flow 
geometry and Reynolds number P,p!0[ . 



Note that the degree of intermittency could be quantified by plotting Dp{r) vs any 
structure function D^{r), (odd order moments taken with the modulus, D*^{r)). 

In this paper we would like to demonstrate that - beyond the mere practical reason of 
better statistics - it seems really essential for physical reasons to take D^^{r) rather than 
D^{r) to have ESS. We do so by examining ESS both for a full numerical simulation |Tl 
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and for Batchelor's parametrization of the structure function |12]. This parametrization will 
give us the opportunity to study finite Reynolds number Re effects. We will furthermore 
discuss the difference in scaling of the longitudinal as compared to the transversal structure 
functions. 

We numerically solve the 3D incompressible Navier-Stokes equation on a A^^ grid with 
periodic boundary conditions. A pseudospectral code is used, the flow is forced on the largest 
length scales. For = 96 we achieve a Taylor- Reynolds number Re\ = 110; scales down 
to 3?7 are resolved; rf = z/^/^/e^^^ is the Kolmogorov length, u the kinematic viscosity, e the 
energy dissipation rate. Time integrations up to 150 large eddy turnovers are performed; 
the flow is locally isotropic to a high degree. More details on the numerical flow are given 
in ref . [Tl[ . 
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FIG. 1. Batchelor parametrization (^) (solid line) of the N=96 data for the longitudinal second 
order structure function (^) (numerical data: circles). We chose a = 12 Arj and L = WSrj. Then 
we used eq. (P) to calculate Djir) (dashed line) which poorly compares with the numerical data 
(squares) for r beyond the VSR. 



II. BATCHELOR'S PARAMETRIZATION 



The longitudinal second order structure function -Df is shown in figure |T]. As in the 
whole paper, lengths are given in multiples of r] and velocities in multiples of {er]Y^^. An 
ISR scaling range is hardly developed because Rex = 110 is still small. The data are very 



well fitted by a parametrization of Batchelor's type |T^-|T6[ with an additional large scale 
cutoff L [|T3|,|r3 
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Here, C = C2 is the asymptotic ISR scaling exponents which from our |11] and others' [0- 
^,||,|18[ ESS analysis we take to be C = 0.70. The only other free parameters are the viscous 
subrange - inertial subrange (VSR-ISR) crossover scale a and, of course, the large scale cutoff 
L. From a fit of eq. (|^) to the present numerical data we find a = 12Ari and L = lOSt]. 

Note that we do not want to imply that all flow fields show a large scale saturation of 
type (^). However, those data (both numerical and experimental) we analyzed (see also 
[|1^,|T^) were well described by eq. (^. As eq. guarantees an analytic behavior of both 
the correction term to on the small scale side and the correction term to r = const on 
the large scale side we think that this is not accidental. 
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FIG. 2. Third order structure function from our numerical simulation (circles), compared with 
the one following from the Batchelor parametrization of -D2 ('') i^^^^ line). The Taylor- Reynolds 
number is Rex = 110. 



With the outer length scale L we can define a Reynolds number Re = ui^rmsL/v. For 
very large Re the structure function (H) develops an ISR scaling law 



l) 



(4) 



Here, 6(2 = C2 ~ 2/3 is the scaling correction to classical scaling and is often called the 
Kolmogorov constant. From our fit we have = e^^^L^'^^a'^-'^^ / {15u) = 2.05. This value 
well agrees with the data b^ = 1.6 — 2.5 known from the literature P-|6|, p!9|] ; Sreenivasan 

ig gives b^ = 2.0 ± 0.4. Note that the fuH structure function D2ir) = 3D^{r) + 
(for isotropic flow) asymptotically scales with a law of type (P, too; the prefactor is 6 ~ 
116^/3 = 6-9. Alternatively, also D2 (r) can be fitted by a Batchelor parametrization 
|13| , |15| with similar quality |15| . 



Our motivation to employ Batchelor's parametrization (^ is to be able to upscale the 
second order structure function -Dg (r) to much larger Re (assuming that b'" and C are fixed 
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at = 2.0 and ( = 0.70) and thereby get consistent data for the transversal second order 
structure function 



D^ir) = {[iuix + r)-uix)).ejr), (5) 

p = 2, and for the third order longitudinal structure function D^{r), which for isotropic, ho- 
mogeneous, incompressible turbulence both follow from (r), namely through the relation 

D^{r)=D^{r) + '-^D^{r) (6) 

and through eq. (^), respectively. 

To upscale -D2 (r) in eq. @, we must know how the parameters a and L depend on the 
Reynolds number. If one accepts Sreenivasan's observations that neither the (asymptotic) 
dimensionless energy dissipation rate [p0|-p2[| = eL/u^.^^^, nor the Kolmogorov constant 

[HH (but note also ref. |2^) depend on the Reynolds number, one gets a weak dependence 
of the VSR-ISR crossover on the Reynolds number |^ , 



Once b^, and Re are fixed, L and Re\ can easily be obtained from above equations |21, 22 
With D^{oo) = 2^2 we get 



.3/8 

) Re"\ (8) 



^ v'^K.^. = ^15 I ] Vi?e. (9) 




Vice versa, once a and L are known, we obtain Re and Rex. For above values the result 



is Re = 520 and Rex = 90 for our numerical flow, in reasonable agreement to the direct 
numerical result Rex = 110. As we will see, the reason for the (modest) underestimation is 
that in the numerical flow there are correlations left at the largest length scales. 



III. dJ and resulting from BATCHELOR'S parametrization 

For -Dj(r) we find poor agreement between the curve evaluated from equation (^) and 
the numerically obtained values, see figure The reason is that at r ~ L there still is 
considerable correlation {u{x + r)u{x)) 7^ 0. More precisely, at the maximal meaningful 
distance r^ax when employing periodic boundary conditions, namely, when r equals half of 
the periodicity length (here, r^ax ~ 146?]), we find {uj{x + e^rmax)uj{x)) / {u'j{x)) ^ 0.25 
and {uj{x+e'^rmax)uj{x)) / {uj{x)) ^ — 0.15 for j = 1,2, or 3. Therefore, D2(^max) is smaller 
and D2{rmax) is larger than 2(m^), which is the value the structure functions would take for 
perfect decorrelation between x and x + Vmax- Geometrically, the above correlations mean 
that there is an eddy with diameter r ~ rmax ~ ttL in the numerical flow. The possibility 
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of such large eddies is a consequence of the periodic boundary conditions (in contrast to 
boundary conditions which put the velocity to zero at the edge of the flow volume) and will 
also survive for larger Rex. Such a large scale eddy implies that the flow is not isotropic and 
homogeneous at the large scales and therefore it should be no surprise that equation 
whose derivation requires isotropy and homogeneity, does not lead to good agreement with 
the data at large scales. Note that the same problem occurs in experimental flow, see e.g. 
figure la of ref. |2^. At the largest measured distance r ~ GSOr; (for Rex = 300) it clearly 
is -Dj(r) > -D2 (r), revealing a large scale eddy. 

Also for D^{r) the agreement between the direct numerical values and those from the 
analytic equation (H), which assumes isotropy and homogeneity, is poor, see figure |. The 
numerical structure function -D3 (r) bends down for large r as the velocity differences at 
large scales have Gaussian hke statistics and consequently odd order moments almost vanish. 
This feature which is not described by eq. (0) is due to the boundary effects; more precisely, 
because there are no larger eddies than on the scale Vmax which could provide correlations. 
In principle, this deficiency can be cured by adding a corresponding term to that equation 
as e.g. done in ref. p5| . 




FIG. 3. Compensated ESS type plot D2/ {D^^)'^/'^ vs D3 for the structure functions from the 
numerical simulation (circles: longitudinal; squares: transversal) and those following from the 
Batchelor parametrization ^ of D2 for Reynolds numbers Re = 5.2 • 10^ as in the numerical 
simulation, for Re = 5.2 • 10^, and for Re = 5.2 • 10^, corresponding to Taylor- Reynolds numbers of 
Rex = 90 (dot-dashed). Rex = 900 (dashed), and Rex = 9000 (solid); C = 0.7, = 2.0. The three 
lower curves are for the longitudinal structure functions, the three upper ones for the transversal 
ones. The arrows indicate lOr/. The external length scales L are beyond the shown regimes. 
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FIG. 4. Compensated ESS type plots £'2/(^3'^)^/^ vs ZJg-^ (three different data sets for three 
different space directions which, however, agree very weU) and D2 / {D'^)'^/^ vs (two data sets, 
also agreeing) for the numerical turbulence, Re\ = 110. The arrows again indicate 10?7; the data 
points are for r = 3r/, r = Gr/, r = 9r/, . . ., left to right. 

IV. HOW TO APPLY ESS? 

We now plot the second vs the third order structure function in a compensated ESS type 
plot, i.e., D^/{D^y/^ vs D^, see figure & The scaling regime in the numerical simulation 
is by far too short to identify any scaling exponent. However, if we repeat this plot, but 
now with D^^ rather than , ESS is seen, see figure ^. The reason is that D2 and D"^^ 
have the same type of large scale saturation (i.e., becoming constant), whereas has a 
different type of large scale behavior (namely, dropping to zero). We have to conclude that 
the extension of the scaling regime by using ESS is mainly an extension towards large scales. 
This even holds if Rex is so small that there is no ISR yet. The small scale onset of scaling 
still is around r ~ lOrj, whether plotting D2 vs D^^ or vs r, which is roughly the crossover 
scale a ~ 12.4// found from employing eq. (|^). In their numerical simulation Briscolini et al. 
find an ESS extension down to r ~ 7?7 (see figure 4 of their paper), roughly the same as 
the 10?7 reported here, but slightly smaller. As pointed out to us by R. Benzi, the origin of 
the slight difference may be that the small scale resolution in Brisolini et al. is down to 
li], whereas here we only have a St] resolution and the deviations in the structure functions 
due to the lower end of the resolution may influence the scaling exponents in a certain range 
of larger scales. 

Our next point is to advocate compensatedESS plots for the visualization of intermittency 
effects. Already Meneveau [0 - see figure 1 of that paper - reveals how misleading an ESS 
plot vs can be. Here, we demonstrate this in figure ^ which shows the original ESS 
plot D2 vs D3, only pretending better scaling than in figure ^. The reason is that in the VSR 
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the scaling exponent is 2/3 for trivial reasons, a value which can hardly be distinguished 
by eye from the ISR value 2/3 + 5^2 ~ 0.70. To avoid this similarity of the VSR and ISR 
exponents, we prefer to use compensated ESS plots |lTTl,p6| . 
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FIG. 5. (a) ESS plots D2{r) vs D^{r) for the longitudinal and transversal structure functions; 
Re\ = 110. The data are the same as in the previous figure, where we showed the compensated 
ESS type plots, in which the different behavior in the VSR and the ISR is clearly visible, 
(b) ESS plots for —D^[r) vs —D^[r). In this figure the Taylor-Reynolds number is only 70, but 
we checked that the lack of ESS does not decrease with increasing Rex. 

If one plots local slopes |^ as done in figure |^a it of course does not make any difference 
whether one takes them from compensated ESS plots or standard ESS plots. From figure §a 
one notices that the ESS scaling ^ 0.69 and C,J ~ 0.72 begins around lOrj. From figure 
lib one also notices that without ESS one could not deduce any scaling exponent at all for 
the small Reynolds number of our numerical calculation. ESS is thus useful already for the 
simple reason that a transition from a local slope of 2/3 to roughly 0.70 is shorter than from 
a local slope 2 to roughly 0.70. 

As we will show now, there is no extended scaling regime towards scales much smaller 
than order ~ IOt], either, if one does ESS type plots with instead of D^^. We do so by 
plotting D2 / {D^y^^ vs with following (via eq. (0)) from the Batchelor parametriza- 
tion eq. (|]) of D2 for various Re, see figure |^. We observe three regimes: The VSR without 
any scaling corrections (i.e., a horizontal line in figure 0), a crossover regime, corresponding 
to the range from r ~ Ir/ to r ~ lOrj, and only for large scales and large Re > 500 the ISR 
scaling corrections 6(2 = 0.033 can be identified. 

To summarize this subsection: There seems to be ESS towards large scales, if the struc- 
ture functions plotted against each other are both calculated with the moduli, i.e., have the 
same large scale saturation behavior. In particular, for the third order longitudinal structure 
function this means that it is essential to take D^^ rather than and to clearly distinguish 
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FIG. 6. (a) Local slopes d Ig /d Ig Df (solid) and d Ig /d Ig L>^'^ (dashed) of the curves m 
figure |5|a. The arrows indicate lOr/. 

(b) Local slopes dlgUJ^ / d\gr (solid) and d Ig L'J/d Ig r (dashed) of the numerical data curves 
D2 ' (r) in figure [l|. 

between the (p and C,p exponents. This was aheady stressed by Stolovitzky and Sreenivasan 
see in particular their figures 2 and 4 where they compare local slopes of vs. and 

^3 ■ 

The natural question to ask is: Is there also ESS for odd order structure functions 
(calculated without the modulus) plotted against each other? In figure [pj we plot — -D5 vs 
—D^. No ESS towards large scales is seen. It seems that odd and even order moments obey 
fundamentally different types of statistics. This finding may be connected to Herweijer and 
van der Water's finding that (p for odd p (calculated from D^{r) oc r^^) are smaller 
than expected from an extrapolation of the neighboring (p^-i (for which the p±l are even). 
A difference between the (p for odd and even p was also found by Stolovitzky et al. P,pm; 
however, for the flow analyzed in that references the (p for odd p are larger than expected 
from the extrapolation. 

We do not understand why nonuniversal forcing and large scale boundary effects roughly 
cancel out in even order structure functions but not in odd order ones (calculated without the 
modulus). With a more elaborate technique it may even be possible to extract nonuniversal 
properties also from ESS plots of even structure functions. On the other hand, we cannot 
exclude that there is more universality in decaying turbulence where less anisotropy through 
the forcing and the boundaries is felt. 
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V. LONGITUDINAL VS TRANSVERSAL STRUCTURE FUNCTIONS 



Next, we focus on the difference in the scahng between longitudinal and transversal 
structure functions, D^{r) oc r^p and Dp{r) oc r'^p , respectively. Recently, different de- 
grees of intermittency for longitudinal and transversal fluctuations were reported in some 
experiments [24, 28, 29| and numerical simulations on decaying turbulence ^0\. We con- 



firmed these findings for statistically stationary turbulence |]Tl| (see also ref. [0]). More 
precisely, it were the ESS type scaling exponents and ^J, defined by oc {D^^)^p and 

oc {Df)^p, which are clearly different; we found 6^^ = 0.21 ±0.01 and = 0.43 ±0.01 
for the deviations from the mean field value ^6 = 2 | ]1 1| | . 

One would be tempted to conclude that the best way to see a deviation in scaling between 
Dp and would be to plot the ratio /Dp vs r (or vs -Dg ). According to eq. (^), -Dj(r) 
and D2{r) scale the same in the ISR, i.e., the ratio should be constant. However, figure ^a 
seems to imply different scaling of -D2 (r) and Dj(r). 

The reason for this apparent discrepancy is that the argument of equal scaling of -D2 (r) 
and Djir) is only valid if both structure functions scale individually. This is not the case 
in the transition ranges or if there is not any ISR yet. Here, the Reynolds number achieved 
in the full simulation is by far too small to give the asymptotic (ISR) scaling. In figure ^ 
we redo this type of plot, but now within the Batchelor parametrization for which we can 
achieve arbitrarily large Re\. Only if Rex > 500 a plateau starts to develop, showing the 
onset of the asymptotically correct ISR behavior. To reliably determine scaling exponents 
from the plateau, one would need at least Taylor-Reynolds numbers ~ 1000 and beyond. 
For Re\ ~ 100 there is a fake scaling law with an apparent exponent of —0.14, which has 
nothing to do with inertial range scaling. 

Going back to eqs. (H) and @, this behavior can be understood. In the VSR we must 
have D'^/D^ = 2 (because of eq. (D and D^ oc r^) and in the ISR we have D^/D^ A/2> 
(because of eq. (H) and roughly D2 oc r^/^), just as seen in figure |^. The crossover between 
these two regimes is about a decade. The same can be seen from figure |^ where besides 
D2 / (D^y^^ we also plotted D2 / {D^Y^^ vs D^. In the crossover regime where the former 
curve bends up, the latter bends down. Again, only for Re\ > 500 the asymptotic scaling 
exponent 5C,2 = — 0.033 starts to be observable. 

The same finite Rex effects which we discussed for the 2nd order structure functions, 
where D2 and are known to have the same scaling, will hinder to determine scaling 
exponents vs r (or vs D^) in higher order structure functions for too low Rex ^ 500. 



VI. SUMMARY 

To conclude, we confirmed the finding of Briscolini et al. 0] that ESS does not extend 
to scales below order ~ 10?]. We furthermore showed from calculations with the Batchelor 
parametrization that scaling exponents C,p calculated from structure functions plotted vs r 
(or vs D^) can only securely be measured for Rex sufficiently larger than 500. For smaller 
Rexi in particular for all present day numerical simulations, one is restricted to relative, ESS 
type scaling exponents calculated from ESS type plots D*^{r) vs D*^{r) and D*^{r) vs 
D*'^{r), whereby it is essential to calculate the structure functions from the moduli of the 
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FIG. 7. (a) The ratio D2 (r) / D2 (r) vs r for the numerical data of a Re\ = 110 (circles) and a 
i?eA = 70 (squares) simulation. Erroneously, one may deduce different scaling oi D2{r) and D2{r) 
because the ratio depends on r. We also included an apparent slope of —0.14. 
(b) D2 (r) / D2 (r) as a function of r for the Batchelor parametrization (^) for Taylor-Reynolds num- 
bers of Rex = 90, Rex = 900, and Rex = 9000. The plateau in the ISR (at roughly lg{A/3) = 0.125) 
expected for isotropic homogeneous flow only starts to develop for as large Rex as Rex ~ 500. The 
VSR plateau is at Ig 2 = 0.30. The r-dependent intermediate range characterizes the transition 
between the VSR and the ISR. The slope is related to its width and the different heights of the two 
plateaus. As a side remark we mention that in this plot one can also notice the Rex dependence 
of the VSR-ISR crossover a, cf. eq. (0): the transition range is shifted towards smaller r with 
increasing Rex- 
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velocity differences. For odd order moments, calculated without taking the modulus, ESS 
does not hold in the presented numerical simulation. 
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